Astronomy & Astrophysics manuscript no. paper 


© ESO 2008 


June 25, 2008 





The infrared emission spectra of compositionally inhomogeneous 
aggregates composed of irregularly shaped constituents 

M. Min 1 , J. W. Hovenier 1 , L. B. F. M. Waters 1 - 2 , and A. de Koter 1 

1 Astronomical institute Anton Pannekoek, University of Amsterdam, Kruislaan 403, 1098 SJ Amsterdam, The Netherlands 

2 Instituut voor Sterrenkunde, Katholieke Universiteit Leuven, Celestijnenlaan 200B, 3001 Heverlee, Belgium 

Received September 15, 1996; accepted March 16, 1997 

Abstract 

Aims. In order to deduce properties of dust in astrophysical environments where dust growth through aggregation is important, knowledge 
of the way aggregated particles interact with radiation, and what information is encoded in the thermal radiation they emit, is needed. The 
emission characteristics are determined by the size and structure of the aggregate and the composition and shape of the constituents. We thus 
aim at performing computations of compositionally inhomogeneous aggregates composed of irregularly shaped constituents. In addition we 
aim at developing an empirical recipe to compute the optical properties of such aggregates in a fast and accurate manner. 
Methods. We performed CDA computations for aggregates of irregularly shaped particles with various compositions. The constituents of 
the aggregate are assumed to be in the Rayleigh regime (i.e. they are much smaller than the wavelength of radiation inside and outside the 
constituents), and in addition we assume that the dominant interaction of the aggregate constituents is through dipole-dipole interactions. We 
computed the spectral structure of the emission efficiency in the lOyum region for aggregates with 30% amorphous carbon and 70% silicates 
by volume with various fractions of crystalline and amorphous components. 

Results. We find that the spectral appearance of the various components of the aggregate are very different and depend on their abundances. 
Most notably, materials that have a very low abundance appear spectroscopically as if they were in very small grains, while more abundant 
materials appear, spectroscopically to reside in larger grains. We construct a fast empirical approximate method, based on the idea of an 
effective medium approximation, to construct the spectra for these aggregates which almost perfectly reproduces the more exact computations. 
This new method is fast enough to be easily implemented in fitting procedures trying to deduce the dust characteristics from astronomical 
observations. 



1. Introduction 

Dust grains in many astronomical environments are assumed 
to be aggregates of small particles of different sizes and com- 
positions. This is especially true for cometary dust and dust 
in disks around young stars. The dust grains in these environ- 
ments are formed by coagulation of smaller particles, leading 
to aggregated structures. In order to interpret observations of 
light scattered, emitted or absorbed by these composite parti- 
cles, a method is required that takes the detailed properties of 
the particles into account. 

Spectral resonances in the refractive index of the parti- 
cle material translate into resonances in the absorption (and 
emission) spectra. Besides composition, the spectral location of 
such resonances is also influenced by the shape and size of the 
particles. Spectra computed for homogeneous spherical parti- 
cles may often display resonances with strengths and spectral 
locations not found in the spectra of irregularly shaped particles 
(see e.g. Fabian et al. 2001 Min et al. 2003) or astronomi- 



cal observations (see e.g. Hony et al. 2002; Bouwman et al. 



2003 Moreno et al. 2003). Much progress has been made 



on computational methods to calculate the optical properties 
of aggregates of spheres (see e.g. Mackowski & Mishchenko 



19961. However, aggregates of homogeneous spherical parti- 



cles are still likely to display resonances typical for homoge- 
neous spheres, and thus behave poorly when comparing to ob- 
servations of natural particles. Test computations we performed 
confirm this. Therefore, a method is required for computing 
the optical properties of aggregates of irregularly shaped con- 
stituents. So far, an efficient and widely employable computa- 
tional technique is still lacking. 

The aggregates expected to be present in astronomical envi- 
ronments like cometary comae or protoplanetary disks are gen- 
erally assumed to be composites of various materials. Although 
the inhomogeneous nature of the particles can be taken into ac- 



count in several ways (see e.g.|Moreno et al. 2003 Schegerer 



|et al. 2006 Kolokolo va et al.j |2007[ ), such inhomogeneous 
structures are usually modeled by assuming separated homoge- 
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2006 ; Min et al.|[2005]|van Boekel et al.||20051[Bouwman et a! 



2003 Brucato et al. 1999| l. Although widely used, the validity 



of this approach is not obvious. In addition it is difficult to in- 
terpret the sizes and shapes derived for the dust grains in this 
way. Therefore, in this paper we study the effects of aggrega- 
tion of various materials in single aggregates on the spectral 
appearance of these particles. In particular we will look at how 
the pronounced spectral structure of crystalline silicates is dis- 
played when they are embedded in a larger aggregate. 

In section [2] we outline the computational approach we 
employ to compute the absorption spectra of inhomogeneous 
aggregates of irregularly shaped particles. Results of calcula- 
tions are presented in section [3] An efficient empirical method 
to compute the spectra for these aggregates using an effective 
medium approach is presented in section|4] Finally, in section|5] 
we present the conclusions of our study. 

2. Computational approach 

Particles very small compared to the wavelength both inside 
and outside the particle, i.e. in the Rayleigh domain, respond 
to electromagnetic radiation as a single dipole with dipole mo- 
ment 

p = aE inc , (1) 

where a is the polarizability of the particle and E; nc is the in- 
coming electric field. The radiation scattered by such a small 
particle is a dipole field, and the mass absorption coefficient, 
i.e. the absorption cross section per unit mass, is given by 



pV 



Im(o'). 



(2) 



Here V is the material volume of the particle, p is the density 
and k - 2n/A where A is the wavelength of incident radiation. 

When we do not have a single small particle, but an ag- 
gregate of such particles, the interaction of all these particles 
with each other has to be taken into account. In other words, 
the incident field each particle experiences is the incoming 
field and the field emitted by all other particles in the aggre- 
gate. In principle all multipole interactions between the con- 
stituents have to be taken into account. This is what is done in 
methods for computing the optical properties of aggregates of 
spheres (see e.g. |Mackowski & Mishchenko] |1996) . Although 
this is in principle not impossible, for irregularly shaped con- 
stituents the computational demand of such an approach is 
very high, and at the moment not feasible. If the constituents 
themselves are very small while most of the interactions take 
place over relatively large distances, like in a large aggre- 
gate of many very small constituents, the dominant interac- 
tion is dipole-dipole interaction. The optical properties of such 
an aggregate of coupled dipoles can be computed using the 
Coupled Dipole Approximation (CD A). In this approximation 
each aggregate constituent is considered to interact as a single 
dipole with a given polarizability. The difference with the well- 



known Discrete Dipole Approximation (DDA, see e.g. Purcell 



& Pennypacker) [1973] |Draine| [1988) is that in DDA the vol- 
ume of the constituents themselves have to be discretized using 
multiple dipoles. This has to be done in order to account for 



the shape of the aggregate constituents. In the CDA we account 
for the shape of the constituents by using for each constituent a 
polarizability representative of an irregularly shaped grain. The 
mathematical formulation of both the CDA and the DDA is ex- 
actly the same; it describes the interaction of dipoles. However, 
the numerical implementation might be slightly different. In the 
DDA one can choose the dipoles to be located on a rectangular 
grid. This allows for the use of Fast Fourier Transform (FFT) 
methods to increase the numerical performance significantly 
( [Goodman et al.| |1991[ jHoekstra et al.| [1998}. Since the con- 
stituents of an aggregate are usually not located on a rectangu- 
lar grid, this is in general not possible when using the CDA. In 
addition, the speedup in DDA using the FFT method requires 
one to fill a rectangular box circumscribing the aggregate with 
dipoles. Since we wish to consider very fluffy grains with a 
low volume filling factor, this is very inefficient and would re- 
quire many additional dipoles. A method providing a similar 
speedup as the FFT method is being developed for the CDA 



computations using the Fast Multipole Method ( Koc & Chew 
[200T]|Amini & Profit||2003| >. 

As we mentioned above the approach of only considering 
the dipole-dipole interactions is valid if the constituents are 
small and the long range interactions between them dominate. 
This is the situation we are currently interested in. Future com- 
putations, using a method taking into account the multipole in- 
teractions, like the superposition T-matrix method or large scale 
DDA simulations using multiple dipoles per constituent, can be 
employed to check the accuracy of the approximation we em- 
ploy here. We wish to consider aggregates consisting of up to 
~ 4000 constituents. In the DDA formulation each dipole has to 
be discretized by, at least, 100 dipoles. Furthermore, the filling 
factor of our aggregates is very low, ~ 5%. This taken together 
implies that for DDA we need ~ 10 7 dipoles to compute the 
aggregates. Although these computations are in principle pos- 
sible on large supercomputers when using an efficient parallel 
DDA program, a single wavelength point will still take up to 
several days computer time. It is therefore not feasible at this 
moment to check the accuracy of the CDA method using DDA 
computations. 

When we consider for each particle of the aggregate only 



dipole interactions we have (Draine 1988 1 



E 



incj 



N 



(3) 



where Pj is the local dipole moment at the position of particle 
j, cij is the polarizability of particle j, Ei nc . ; is the incoming 
field at the location of particle j, and N is the total number of 
constituents in the aggregate. The matrix A;* determines the 
electric field at the position of particle j due to the dipole field 
emitted by particle k (for details see e.g. Draine |1988|l. When 



Eq. ([3]) is solved for the P ; the mass absorption coefficient can 
be computed by 



Ank 



pV\E n 



N ( 2 1 

J]|lm[P r (a 7 1 rP*]--^P}.P j |, (4) 
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where V is the total material volume of the aggregate, p is the 
average material density, and the asterisks denote the complex 
conjugates. 

The dipole polarizability of a homogeneous spherical par- 
ticle is given by 



3V- 



m 2 - 1 



+ 2 



(5) 



where m is the complex refractive index of the particle material. 
The polarizability of irregularly shaped particles depends on 
the particle shape. However, a reasonable approximation for 
the average polarizability of an ensemble of highly irregular 
particles is given by the Continuous Distribution of Ellipsoids 
(CDE; see |Bohren & Huffman! [1983) |Min et aL||2006b| >. The 
CDE polarizability is given by 



a = 2V 



m 2 — 1 



In m" 



1 



(6) 



More realistic polarizabilities can be computed using the 
method of |Min et aL] ( |2006b| l. However, here we wish to 
study only the effects of using irregularly shaped particles as 
monomers so we employ the frequently used CDE polarizabil- 
ity for the aggregate constituents. 

In the following subsections we will first define two quanti- 
ties we use for the analysis of the emission spectra of the aggre- 
gates. Then we will discuss how the aggregates are constructed 
and what composition we will assume. 

2.1. Effective mass absorption coefficients 

The mass absorption coefficient as given by Eq. Q is basically 
the sum over the absorption contributions of each component 
of the aggregate. We can thus easily compute the effective mass 
absorption coefficient of a part of the aggregate by summing 
over only that part of the components. In this way we can cre- 
ate an effective mass absorption coefficient for each material 
component. The effective mass absorption coefficient of mate- 
rial component i of the aggregate is defined as 



ftpiV\E h 



^ jlm [P; • (a^yP*] - ^ 3 P* • Pjl , (7) 



where is the volume fraction of the total aggregate com- 
posed of material component i, pj is the material density of this 
component, and the summation over j now only involves the 
dipoles composed of material i. When there are N mat different 
materials that make up the aggregate, the total mass absorption 
coefficient is simply given by 



K = > Ki . 

P 



(8) 



2.2. Construction of the aggregates 

Aggregates in astronomical environments form by coagula- 
tion of small grains. The structures that form in this way can 
have various shapes and degrees of compactness depending on 
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Figure 1. The real (solid) and imaginary (dotted) part of the 
polarizability per unit volume for the various components in- 
cluded in the aggregates in the 10 pun region. For the forsterite 
and enstatite we display the polarizability per unit volume av- 
eraged over the three crystallographic axes. 



the environmental conditions. On the one hand, fluffy aggre- 
gates might form under conditions where aggregates of approx- 
imately equal sizes coagulate to form larger structures (Kempf 



et al. 1999[ ). On the other hand, very compact structures might 



form when single monomers are added to a larger aggregate 
( |Ball & Witten] [1984) . As shown by |Min et al.| ( |2006a| l the 
compactness of the aggregate has a large impact on its spec- 
tral appearance. For very fluffy aggregates the spectral struc- 
ture typical for their constituents is visible, whereas for very 
compact aggregates the spectral structure of the total aggregate 
dominates. 

We construct the aggregates by adding single particles to 
a growing aggregate. We start out with a single particle and 
shoot a second particle with a random impact parameter at the 
particle till it hits the particle and sticks. The resulting cluster 
is randomly rotated and a new particle is shot at the aggregate 
in the same way. In this way we produce moderately fluffy ag- 
gregates. In order to construct aggregates of irregularly shaped 
particles we shoot Gaussian Random Field (GRF) particles at 



each other (Grynko & Shkuratov 2003 Shkuratov & Grynko 



2005]), thus constructing an aggregate of these very irregularly 
shaped particles. We then locate a dipole with CDE polariz- 
ability at the center of mass of each GRF particle. The volume 
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Name 


Composition 


Lattice Structure 


Volume fraction 


Ref. 


Amorphous silicate mixture 


Fe rich Olivine 


MgFeSi0 4 


Amorphous 


8.6% 


[1] 


Mg rich Olivine 


Mg 2 Si0 4 


Amorphous 


47.7% 


[2] 


Mg rich Pyroxene 


MgSiOj 


Amorphous 


37.0% 


[1J 


Na/Al Pyroxene 


NaAlSi 2 6 


Amorphous 


6.7% 


[3] 


Other components 


Forsterite 


Mg 2 Si0 4 


Crystalline 




[4] 


Enstatite 


MgSi0 3 


Crystalline 




[5] 


Carbon 


C 


Amorphous 




[6] 



Table 1. The materials used in our aggregates. In the case of the amorphous silicates we use the names 'olivine' or 'pyroxene' 
to indicate the average silicate stoichiometry. The volume fractions of the amorphous silicates used to construct the average 
amorphous silicate polarizability are also listed (taken from Min et al. 2007| l. Note that these volume fractions represent the 



relative contributions of these species to the total silicate component in the aggregate. Also note that these fractions are given as 
volume fractions rather than mass fractions as in |Min et"aL] ( |2007] l. The references in the last column refer to: [1] Dor schner et al. 



(fT995ll, [2]|Fabian et al.| f2000), [3] |Mutschke eFaT1 ( fT998) , [4] |Servoin & Piriou| ( [T973] >, [5] |Jageret al^ ( [T998) , and [6] |Preibisch 
le t al.| ( |1993| >. ~ 



equivalent radius of the aggregate constituents is taken to be 
0.4/im. 



2.3. Composition 

We consider inhomogeneous aggregates. We take 30% of the 
constituents of the aggregate to be composed of amorphous 
carbon and 70% to be composed of silicate. All constituents 
have the same size, so these, and any further abundances men- 
tioned in the text, are volume fractions. When we assume that 
all of the Si-atoms are in the silicates, the carbon/silicate ratio 
we take corresponds to approximately half of the solar abun- 
dance of carbon atoms to be in the grains. Of the 70% silicate 
a fraction is crystalline while the rest is amorphous. Crystalline 
silicates display a very pronounced spectral structure, and part 
of our study is to see how this spectral structure is displayed 
for various fractions of crystalline silicates. 

The polarizabilities are computed using Eq. (|6]l with labo- 
ratory measurements for the refractive index as a function of 
wavelength. For the amorphous carbon we used measurements 



by Preibisch et al. ( 1993 ). Forsterite and enstatite are crystalline 
materials that have a refractive index different for the electric 
field oriented along the different crystallographic axis. We take 
this into account in the computations coupling the dipoles. The 
refractive index data was taken from |Servoin & Piriou| ( fT973] l 
and |Jager et al. (1998 ) for the forsterite and the enstatite re- 
spectively. 

For the polarizability of the amorphous silicate constituents 
we take a mixture of various amorphous silicates according 
to the composition derived for the interstellar amorphous sil- 
icates by Min et al. ( 2007 1. This mixture is given in table [TJ 
and represents an average amorphous silicate composition of 
Mg ] .36Feo.13Nao.05Alo.05SiO3.59 ■ 

The polarizabilities as functions of the wavelength for all 
different materials used in our computations are shown in 
Fig.[T] 



3. Results 

In Figs. [2] we plot the total mass absorption coefficient for 
the various aggregates. Spectra are shown for various volume 
equivalent radii. By volume equivalent radius we mean the ra- 
dius of a sphere with an equivalent volume as the aggregate. 
We clearly see the expected flattening of the 10 fim feature with 
increasing aggregate size. Also, for the more crystalline aggre- 
gates we see the spectral signature of the crystalline species 
appearing. 

In Figs. [3] and |4] we plot the effective mass absorption co- 
efficients of the amorphous silicate (Fig. [3} and crystalline 
forsterite (Fig.|4]i components. In these plots we also show the 
mass absorption coefficients for the cases in which the entire 
aggregate is composed of a single material (amorphous silicate 
and crystalline forsterite respectively). 

The abundance of amorphous silicate material in the aggre- 
gate is quite large. Therefore, it is not very surprising that for 
the amorphous silicate the pure and effective mass absorption 
spectra are quite similar (see Pig. [3]». For the forsterite compo- 
nent this is different. As can be seen from Fig. [4] the effective 
mass absorption coefficients show a spectral signature typical 
for small grains even when the aggregate as a whole is large. 

From these results we conclude that small inclusions of ma- 
terials with strong resonances can show the spectral signature 
of small grains even when they are embedded in a larger grain. 
This is an important result to keep in mind when interpreting 
infrared spectra of environments where coagulated grains are 
expected to be present. 

To quantify the statement above we determined the appar- 
ent size for the various components making up the aggregate. 
This is the size that the components appear to have spectro- 
scopically. To do this we applied the following procedure. 

1. We computed the effective mass absorption coefficients, 
of the various materials as function of wavelength. 

2. We computed the mass absorption coefficients for homoge- 
neous aggregates, i.e. consisting of a single material, for all 
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Figure 2. The total mass aborption coefficients for the aggregates with crystalline silicate fraction (forsterite + enstatite) of 2% 
(left panels), 10% (middle panels), and 20% (right panels). The volume equivalent radii of the aggregates are 1 /mi (upper panels), 
3 //m (middle panels), and 7 /mi (lower panels). 
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Figure 3. The effective mass absorption coefficient of the amor- 
phous silicate component in the aggregates (full lines). The 
plots are organized as in Fig. [2] As a reference we also plot 
the mass absorption coefficients for an aggregate with the same 
size but composed of pure amorphous silicate, i.e. no carbon or 
crystalline silicates in the aggregate (dashed line). 
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Figure 4. The effective mass absorption coefficient of the crys- 
talline forsterite component in the aggregates (full lines). The 
plots are organized as in Fig. [2] As a reference we also plot 
the mass absorption coefficients for an aggregate with the same 
size but composed of pure crystalline forsterite (dashed line). 



materials and for all sizes up to 7/<m for which we have 
computations (1, 2, 3, 4, 5, 6 and 7 /mi). 
3. The spectral shape of the effective mass absorption spectra 
computed in step 1 were compared with the mass absorp- 
tion spectra for homogeneous aggregates computed in step 
2 to determine the apparent sizes for all materials in the ag- 



gregate. In this comparison we also allowed a scaling factor 
for the spectra to be fitted. In this way we fit the spectral 
shape and not to the absolute value of the mass absorption 
coefficient. 

The results of this procedure are shown in Fig. [5] It is clear from 
this figure that the forsterite and enstatite grains appear much 
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Amorphous silicate 
Forsterite 
Enstatite 




Inhomogeneous aggregate size (/im) 

Figure 5. The fitted apparent size for the different materials in- 
side the aggregate as functions of the aggregate size for a crys- 
tallinity of 20%. Results for crystallinities of 2 and 10% show 
similar trends. Along the horizontal axis we plot the volume 
equivalent radius of the inhomogeneous aggregates. Along the 
vertical axis we plot the volume equivalent radius of the best fit 
homogeneous aggregate. For example, we can consider an in- 
homogeneous aggregate with the composition as given in the 
text, a crystallinity of 20%, and a volume equivalent radius 
of 4 fim. From this figure we can now determine that this ag- 
gregate is best represented spectroscopically by a polydisperse 
mixture of homogeneous aggregates with amorphous silicate 
aggregates having a volume equivalent radius of ry = 3.3 /mi, 
forsterite aggregates with ry — 1.3 fim, and enstatite aggregates 
with r v = 1.9 //m. 

smaller than the amorphous silicate grains. This is caused by 
the small volume fraction of the crystalline materials. 

Note that although the spectral shape of the low abundance 
materials is representative for small grains, the strength is re- 
duced compared to the mass absorption coefficient of individ- 
ual small grains. This makes it very hard to get reliable abun- 
dance information when using single grains to analyze aggre- 
gate spectra. In the next section we will discuss an approximate 
method that can be used to correctly obtain abundance informa- 
tion from absorption or emission spectra caused by aggregate 
dust particles. 

4. Effective medium approach 

In this section we outline a fast, empirical, effective medium 
approach to simulate the absorption spectra of the aggregates 
discussed above. In contrast to the CDA computations, which 
can take up to several days to obtain a spectrum, this method 
requires only Mie computations which can be done for any 
particle size in negligible computer time. As we will show, 
this method almost perfectly reproduces the CDA computa- 
tions while the computation time required is much smaller. The 
method we propose is based on the idea of simulating the en- 
tire aggregate as a sphere with an effective refractive index. 
There are several methods in the literature to compute the ef- 
fective refractive index under certain assumptions. Recently 



fective medium theories for the optical properties of porous 
particles in the visual part of the spectrum. In our case the ag- 
gregates are very fluffy, open structures. If one would draw a 
circumscribing sphere around the aggregate the vacuum frac- 
tion inside this sphere would be very large. 

The two most commonly applied mixing rules, the Garnett 
and Bruggeman mixing rules, can be derived from the same 
basic equation 



a a {m eS lm m ) = fm\^fia c (mi/m m ) + (1 - /fiii)ar c (l/m m ), 



(9) 



where a a (m) and a c (m) are the polarizabilities of the aggregate 
as a whole and the components of the aggregate respectively, 
m e ff is the effective refractive index of the material, m m is the 
so-called matrix material, /; is the volume fraction of compo- 
nent i, m, is its refractive index, and /mi is the filling factor of 
the aggregate. Note that also vacuum is included as a compo- 
nent, with mi = 1 and volume fraction (1 - fm)- The difference 
between the Garnett and the Bruggeman mixing rule is in the 
choice of the matrix material m m . This refractive index repre- 
sents the material in which all components are assumed to be 
suspended. In the Garnett rule the dominant component is used 
for m m . However, when none of the components is clearly dom- 
inant, a better choice is to take m m = m e g, so all components 
are assumed to be suspended in the effective medium itself. This 



results in the Bruggeman mixing rule. Indeed, Voshchinnikov 



et al. ( 2007| l find that the Bruggeman mixing rule generally 
gives better results. However, when one of the materials clearly 
dominates, both mixing rules converge towards each other. In 
our case the filling factor of the aggregates is so low that vac- 
uum is clearly the dominant component (~95%). Therefore, we 
use the Garnett mixing rule taking m m = 1 . 

In the classical form of the mixing rules both a a and a c 
are taken to be the polarizabilities of spheres. For irregularly 
shaped constituents we can adapt this by taking in Eq. (|9]l the 
polarizability a c to be that of an irregularly shaped particle. The 
aggregate as a whole is still approximately spherical so for a a 
we take the polarizability of a sphere. Also substituting m m = 1 
this results in 



3 V—; = /mi > fi a c (mi). 

m „ + 2 

eft (=1 



(10) 



In our case we use the CDE polarizability for the a c as given 
by Eq. |6]). We note here that one could also have written the 
polarizability for a CDE grain on the left hand side, i.e. for 
a a . Since the aggregates are very fluffy structures m^s is gener- 
ally close to unity, so the effect of the circumscribing shape of 
the aggregates (which is what a a represents) is generally very 
small. Thus, in order to do computations for different aggregate 
sizes one might simply increase the size of the circumscribing 
sphere and perform Mie computations using the effective re- 



fractive index as given by Eq. ( 10 1. We will refer to this method 
as the Aggregate Polarizability Mixing Rule (APMR). 

Now we have to determine how to get the required /mi and 



Voshchinnikov et al. (2007i have analyzed several of these ef- radius of the homogeneous sphere. As is shown by Min et al. 
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Figure 6. A comparison of the spectra computed using the CDA computations and the effective medium approximation (APMR) 
for various values of the crystallinity (forsterite + enstatite volume fraction). The solid line represents the CDA computations 
while the dashed line respresents the resulting opacities from the effective medium approximation. Below each panel the relative 
error between the two curves is shown. Scale in these panels runs from -10% to +10%. The differences are in all cases smaller 
than 6%. 



(2006a) the most natural choice, namely using the circumscrib- 
ing sphere, leads to an overestimate of the result of the fluffy- 
ness of the aggregate. Computations we performed for the case 
at hand confirm this conclusion. This is because the aggregate 
constituents are not randomly distributed in the volume circum- 
scribed by this sphere. Therefore, we have to choose a some- 
what smaller radius which also somehow takes into account the 
structure of the aggregates. A choice, which leads to excellent 
results as shown below, is to use the radius of gyration of the 
aggregate. 



1 w 



r„| 2 



(ID 



In this equation r, is the location of constituent i, and r ( > is the 
center of mass of the aggregate. For fractal aggregates with 
fractal dimension Df the radius of gyration can be expressed 



as (Filippovetal. 2000) 



r g=7 r v > 



(12) 



where y is a constant depending only on the size of the con- 
stituents and on the so-called fractal prefactor. By comparing 



the gyration radii computed using Eqs. (Hi and ( 12 1 one finds 
that the aggregates we use are best represented using Df — 2.82 
and y = 2.44. The filling factor is now defined by 



ff 



fin 



(13) 



where ry and r g are the volume equivalent radius and the gy- 
ration radius respectively. For the radius of the homogeneous 



sphere with the effective refractive index as given by Eq. ( 10 1 
in the final computation we now have to take r g , which results 
in the aggregate and the sphere having the same mass. 

Substituting m m = m e g, the Bruggeman rule for irregularly 
shaped constituents results in 

/V 

/fill ^ fi a c ('«,7'«eff) + (1 - /fiii)a c (l /'«eff) = 0. (14) 
(=1 

The advantage of using the Garnett mixing rule over the 
Bruggeman rule is that it has an analytic solution for m e ff while 
the Bruggeman rule requires one to solve Eq. ( 14 1 to find m e g . 



We have performed computations using the Bruggeman rule 
and find that in our case the differences with the Garnet rule 



according to Eq. ( 10 1 are negligible. As explained above, this is 
because the aggregates we consider generally have a filling fac- 
tor of less than 5%. If one would like to extend the equations for 
the APMR to aggregates of different structures this is not neces- 



sarily the case. Considering the results of Voshchinnikov et al. 
(2007]l it is to be expected that especially for aggregates with a 
larger filling factor, i.e. lower porosity, it might be wise to use 
the Bruggeman mixing rule according to Eq. ( 14 1 in stead of the 



Garnett rule (Eq. 10 1. However, in general aggregates formed 
in protoplanetary disks are expected to have fractal dimensions 
lower than the ones we used and thus be even more porous. 

With the simple equations above we can almost perfectly 
reproduce all our CDA computations. In Fig. [6] we compare 
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the spectra for a mixture of 2, 10 and 20% crystalline silicates 
computed using the CDA method outlined in section[2]with the 
effective medium computations as outlined above. As one can 
see, the match is almost perfect. 

5. Conclusions 

We have studied the spectral behavior of inhomogeneous ag- 
gregates of irregularly shaped constituents. We find that the 
spectral appearance of such aggregates is a complex interplay 
between the abundances of the materials and the size of the ag- 
gregate. Also the shape and structure of the aggregate plays an 
important part in the spectral appearance of such particles. 

The most striking conclusion of our computations is that 
the apparent particle size of a dust component, i.e. the size 
derived from the infrared spectrum, is a strong function of its 
abundance. Materials with a relatively low abundance will ap- 
pear spectroscopic ally smaller than materials with a relatively 
high abundance. 

We present a very efficient, fast and accurate approximate 
method to compute the opacities of the inhomogeneous aggre- 
gates using an empirical effective medium approach. We refer 
to this approach as the Aggregate Polarizability Mixing Rule 
(APMR). This approach can be used to study the effects of in- 
homogeneity of particles of larger sizes than can be computed 
using more exact methods given the large computational de- 
mand that these methods put on the computations for increas- 
ing particle size. 

Often, the analysis of infrared spectra of dusty environ- 
ments is performed by assuming homogeneous, solid grains 
with various compositions and a limited number of grain sizes. 
The APMR approximate method is fast enough to be easily 
implemented in fitting algorithms trying to deduce the compo- 
sitions and sizes of dust particles in these astrophysical envi- 
ronments. Furthermore, we propose to use the APMR to gauge 
the results obtained by previous studies and further investigate 
the effects of grain aggregation on the analysis of infrared spec- 
tra. Results of such an analysis will be presented in a separate 
paper. 
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